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ABSTRACT 

We present quantitative morphology measurements of a sample of optically se- 
lected group galaxies at 0.3 < z < 0.55 using the Hubble Space Telescope (HST) Ad- 
vanced Camera for Surveys (ACS) and the gim2d surface brightness-fitting software 
package. The group sample is derived from the Canadian Network for Observational 
Cosmology Field Redshift survey (CN0C2) and follow-up Magellan spectroscopy. We 
compare these measurements to a similarly selected group sample from the Millen- 
nium Galaxy Catalogue (MGC) at 0.05 < z < 0.12. We find that, at both epochs, 
the group and field fractional bulge luminosity (B/T) distributions differ significantly, 
with the dominant difference being a deficit of disk-dominated (B/T < 0.2) galax- 
ies in the group samples. At fixed luminosity, z=0.4 groups have ~ 5.5 ± 2 % fewer 
disk-dominated galaxies than the field, while by z=0.1 this difference has increased to 
~ 19 ± 6 %. Despite the morphological evolution we see no evidence that the group 
environment is actively perturbing or otherwise affecting the entire existing disk popu- 
lation. At both redshifts, the disks of group galaxies have similar scaling relations and 
show similar median asymmetries as the disks of field galaxies. We do find evidence 
that the fraction of highly asymmetric, bulge-dominated galaxies is 6 ± 3 % higher in 
groups than in the field, suggesting there may be enhanced merging in group environ- 
ments. We replicate our group samples at z = 0.4 and z = using the semi-analytic 
galaxy catalogues of IBower et al.l ((2006|). This model accurately reproduces the B/T 
distributions of the group and field at z=0.1. However, the model does not reproduce 
our finding that the deficit of disks in groups has increased significantly since z=0.4. 
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1 INTRODUCTION 

Galaxies, at a simple level, are a mixture of two fundamen- 
tal and distinct components: a bulge and a disk. In the local 
universe, bulge dominated galaxies are generally red and 
quiescent, while disk dominated g alaxies are ge nerally blue 
and actively forming stars (Blanton et al ] 120031 ). Thus, the 
morphology of galaxies may be important when trying to 
explain the observations that show the cosmic star forma- 
tion density has r apidly declined fro r a a peak at z ^l-1.5 
(jLillv et al.l Il996l : iMadau et all Il996l : iHopkinsI [2004 ). and 
that the fractio n of red galaxies has rapidly incre ased over 
the same time (jPaber et al.ll2007l : iBeTl et al.ll2004h . Indeed, 
these observations suggest that the process which transforms 



galaxies from disk-dominated to bulge-dominated is also the 
process which tra nsforms them fr om the blue cloud to the 
red sequence feg. jBell et n]l2007ll. Howev er, studies of red, 
disk-dominated gala xies (IWolf et al.ll2005l) and blue, bul^e- 
dominated galaxies ([Abraham et al.lll991 : iMenanteau et al.l 
2006) suggest this model may be too general. In addition, 
the large fractions of passive spirals at intermediate red- 
shift suggests that the truncation of star formation may 
happ en before the mor phologic al transform ation mechanism 
(Pogg ianti et al.lll999l . but see IWilman~ al. 2008). 

The local environment of a galaxy is an important fac- 
tor in its evolution. Observations of galaxy clusters have 
shown that galaxies within clusters ha ve lower star forma- 
tion rates than the general field (eg. iBalogh et al.l ll999l V 
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The rapid structure growth associated with ACDM cosmol- 
ogy, which increases the local galaxy density of the average 
galaxy with time, may be the key driver of the decline of 
galaxy star formation rates. Indeed, analogous to the lower 
star formation rates in clusters, there is a correlation be- 
tween the local galax y surf ace density and the morphology 
of galaxies (Dressier Il980l ). At low redshift, the percent- 
age of early type galaxies increases, and the percentage of 
late types dec reases, with increasing density. Interestingly, 
iDresslerl (|l980l l found that this relation was equally strong in 
centrally-concentrated, relaxed clusters and in irregular, less 
centrally-con centrated clusters. At higher redshift, z ^ 0.5, 
^Dressie r et al.l (|1997i) showed that this morphology- density 
relation is stronger in highly concentrated clusters than in 
less concentrated clusters. 

Although these studies point to the crucial role clusters 
play in the morphological transformation of galaxies, they 
are rare environments, and thus cannot have a large enough 
effect on the properties of galaxies to explain the decline 
of the star formation density of the universe as a whole. 
However, the less dense environment of optically selected 
groups is the most common env ironment for galaxies in the 
local universe Indeed, IPostman &; Gellerl 

("1984') found that the morphology- density relation extends 
smoothly into the group scale environments. Further sup- 
porting the integral role of groups, suppressed galaxy star 
formation rates in group-scale environme nts of the loca l 
universe is now wel l-established (eg., Balo gh et al.l l20o3 l. 
IWil man et al.' ('2005a) have shown that the fraction of galax- 
ies with [OII]A3727A emission, a measure of star formation, 
is much higher in group galaxies at intermediate redshift, z 
'-^0.4, than in the local universe; however, the group galaxies 
still exhibit suppressed star formation relative to the field at 
the same epoch. 

The physical cause of the suppression of star formation 
since z^-^l isn't clear, but there are many candidates, each 
with their own morphological signatures. Within the con- 
text of ACDM cosmology, gala xy mergers are often thought 
to be a dominant mechanism (|HoDkins et al. I l2007l ). simu- 
lations suggest that a major merger between two gas- rich 
and star-forming spiral galaxies produces a gas-poor, pas- 
sively evolving elliptical galaxy (|Toomre Toomrd Il972l : 
iMihos k HernauistllT996l ). If dominant, this scenario sug- 
gests that quiescent spiral galaxies should be rare, and that 
the transformation of morphological type should precede or 
happen at the same time as the complete suppression of 
star- format ion. Group environments are thought to be the 
ideal place for galaxy mergers because of their high density 
and small relative velocities. 

Recently, driven by dual observations of large bub- 
bles seen in the hot X- r ay gas of the intracluster medium 
(|McNamara et al.1 |2000l: iFabian et all 2000) and the cor- 
relation between the mass of the galacti c bulge and the 
size of the central supermassive black hole (jMagorrian et all 
ll998l : lFerrarese k Merrit"tll2000h . feedback from active galac- 
tic nuclei has become a popular explanation for the sup- 
pression of star formation rates in massive galaxies. Semi- 
analytic galaxy formation models have succe ssfully intro- 
duced these mechanisms in a parametrised way (|Bower et al.l 
I2OO6I : [Croton et al]|2006l l, but the details are still uncertain. 
Such energy feedback mechanisms may not directly alter the 



galaxy morphology, but reduced star formation may result 
in significant fading of the disk component. 

Meaningful morphological measurements are necessary 
to break the degeneracy of physical explanations of star for- 
mation truncation. Visual classification of galaxies onto a 
Hubble (or similar) system has proven to be very useful 
for the study of galaxy evolution. However, the high resolu- 
tion and uniform quality of large galaxy surveys has given 
rise to automated morphology systems which attempt to 
make more quantitative measurements than a visual sys- 
tem will allow. Non - parametric morp hology systems (eg. 
[Abraham eraD l2003l : [Lotz et all 120041 ) are robust, but are 
not easily linked to physical quantities such as bulge or disk 
scale lengths. For this reason, in t his paper we use a popular 
code, gim2d (jSimard et al.|[20o3 ). to fit parametric models 
to the surface brightness profiles. Parametric systems suf- 
fer because they fit an a priori model to the galaxy surface 
brightness and, as such, are prone to giving non-physical 
results in some cases. We therefore adopt the logical filter- 
ing system proposed by I Allen et al.l (|2006l ) , to help mitigate 
some of these effects. 

In this paper, we examine the morphological properties 
of optically selected samples of group galaxies at z=0.4 and 
z=0.1. In 32] we describe our data samples and our morpho- 
logical measurements. In 33]we present the main data results 
and in 311 we discuss what the data tells us about galaxies in 
transformation and compare our data resu l ts wit h the semi- 
analytic galaxy catalogue of iBower et al.l (|2006h . We sum- 
marize our main results in JS] Details of the group-finding 
algorithm are given in Appendix A, while in Appendix B we 
consider possible systematic effects on our results. Finally, 
in Appendix C we show a representative sample of images 
from our CN0C2 z—OA sample of galaxies. Throughout this 
paper we assume a cosmology with matter density =0.3, 
energy density Qa =0.7, and present-day Hubble constant 
Ho= lOO/i km s"^ Mpc"^ with h = 0.75 (or /its = 1). 



2 THE DATA 

2.1 The 0.3 ^ z ^ 0.55 Sample 

Our moderate-redshift galaxy sample is derived from the 
Canadian Network for Observational Cosmology Field 
Galaxy Redshift Survey (CN0C2), a spectroscopic and pho- 
tometric survey completed with the Multi-Object Spectro- 
graph instrument at th e 3.6-m Canada France Hawaii tele- 
scope (|Yee et al.ll2000l l. The survey was designed to study 
galaxy clustering and evolution. It targeted galaxies in the 
redshift range 0.1 < z < 0.6 over four different patches of 
sky totaling about 1.5 square degrees. The survey consists 
of 5 colour (U,B,V,Rc, Ic) photometry of '-^40,000 galax- 
ies to a limiting magnitude of Rc=23.0 mag. Spectroscopic 
redshifts of ^^6000 galaxies were obtained with an overall 
sampling rate of 48% to Rc=21.5. This large survey allowed 
ICarlberg et al. 1 (|200lh to identify a set of 200 groups using 
a friends- of- friends redshift-space group finder. 

Wilman et al.1 (l2005bl ) followed the CN0C2 survey with 
deeper spectroscopy of a set of 26 groups (2 targeted, 6 
serendipitous) drawn from the ICarlberg et al.l catalogue us- 
ing the Low Dispersion Survey Spectrograph (LDSS2) at 
the 6.5m Baade telescope at Las Campanas Observatory in 
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Figure 1. We show k-correcti ons to rest frame V band magni- 
tudes, derived from a suite of iBruzual CharlotI {2003) popu- 
lation synthesis models. The solid lines indicate a constant star 
formation rate; the dashed lines are models with an exponentially 
declining star formation rate; and the dotted lines represent single 
stellar population models created at z = 4. All models are shown 
with and without one magnitude of dust extinction. The thick 
red line indicates the adopted k-correction of = —0.75 2; + 0.65. 



Chile. These groups were chosen to lie within 0.3 <z< 0.55, 
and galaxies brighter than Rc=22.0 were targeted for spec- 
troscopy. This additional spectroscopy was designed to give 
near full c ompleteness at bright magnitudes (Rc<20; for 
details, see IWilman et al.ll2005bh . 

For each of the 20 targeted groups, we obtained single 
orbit Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) pointings in the F775W filter during Cycle 
12. These data were proces s ed wi th the ACS pipeline as de- 
scribed bv lPavlovskv et. all (|2005l l . The images were further 
processed with the Multidrizzle task in pyraf to remove cos- 
mic rays and hot pixels. 

Sources were detected in the HST ACS images using the 
SExtractor software v2.3.2 (Bertin & Arnouts 1996). For a 
source to be accepted, the signal in at least 10 of its ACS 
pixels (0.5 arcsec^) had to be a minimum of 1.3 a above 
the background. The faintest sources that are reliably de- 
tected using these criteria have R775 ~ 23.9. In this paper 
we restrict the analysis to sources with redshifts (Rc < 22) 
and are therefore insensitive to these detection parameters. 
However, we are sensitive to the deblending parameters as 
the automated surface brightness fits use the segmentation 
image produced by SExtractor to identify which pixels be- 
long to the galaxy. We used 32 deblending subthresholds, 
with a minimum contrast parameter of 9.0 x 10~^. By trial 
and error, these parameters gave the best deblending upon 
visual inspection of the output segmentation im ages. 

We have used a suite of iBruzual CharlotI (| 20031 ) pop- 
ulation synthesis models to k-correct our HST magnitudes. 
We have chosen to avoid large k-corrections by correcting 
all galaxies to the nearest restframe waveband, V. Figure [1] 



shows the models used. The k-correction is mainly sensitive 
to the star formation history and dust attenuation of the 
galaxy, and is insensitive to the initial mass function and 
the metallicity. We add a correction of k' = —0.75 2; + 0.65, 
shown by the red dashed line, to each of our galaxies. This 
was chosen to minimize bias for any single galaxy type, but 
dominates the uncertainty in the magnitudes to =b 0.1. In all 
cases, our statistical uncertainties dominate over this source 
of systematic error. All HST magnitudes are quoted in the 
k- corrected rest frame V band. 

2.1.1 The Group and Field Samples 

The group catalogue of ICarlberg et al.l (|200lh originally 
identified virialized galaxy groups in redshift space using 
an iterative friends-of-friends algorithm on the CN0C2 red- 
shift catalog. They found over 200 groups, with an average of 
3.8 confirmed members per group. To take advantage of the 
deeper and more complete LDSS2 spectroscop y available to 
us, we redefine the CN0C2 groups following Wilman et al.l 
(2005b). We discuss the group finding algorithm in more de- 
tail in Appendix |A| Briefly, an iterative procedure is used, 
which initially selects galaxies within two times the velocity 
dispersion of the mean group redshift, and with a transverse 
distance from the group centre within 1/5 of the dispersion 
distance. In each iteration, the velocity di spersion is re com- 
puted using the Capper estimator (Beers et al.lll99Cil l and 
the centre is recomputed as the luminosity- weighted geomet- 
ric centre of the group. The iterations are continued until a 
stable group membership is reached. 

Using these group centers and velocity dispersions, we 
restrict the membership of our group sample to galaxies 
within two velocity dispersions of the group redshift and 
within 500/175^ kpc of the group center in the transverse di- 
rection. To obtain a true sample of group-sized halos we fur- 
ther restrict our group sample to have velocity dispersions 
< 700 km/s within 500h^^kpc of the group center. 

Field samples of galaxies are commonly defined in one of 
two ways: as an "isolated" sample, in which galaxies within 
groups and clusters are removed, or as a "global" sample, 
which includes all galaxies regardless of their environment. 
In practice, the removal of all group and cluster galaxies 
is not possible in our sample. Incomplete redshift sampling 
and the observational uncertainties associated with group 
membership would lead to an "isolated" field sample which 
still contained some group galaxies. Therefore, we prefer to 
define a field sample that contains all galaxies, regardless of 
their group or cluster membership. However, our follow-up 
spectroscopy focused on regions that have groups identified 
in the CN0C2 survey, and our morphologies are derived 
from small ACS images centered on our each of 20 targeted 
groups. To avoid the bias toward groups that would other- 
wise be present, we define our field sample to include only 
those galaxies that are not in groups as identified by the 
Carlberg et al. algorithm. As discussed in ^2.2.11 and Ap- 
pendix [A] because of the incompleteness of the CN0C2 sur- 
vey, and the strict group finding algorithm of Carlberg, this 
leaves a field sample which is only slightly depleted in group 
galaxies when compared to the universe as a whole at z=0.4. 
A similar selection in the semi-analytic group catalogues dis- 
cussed later f ^4.3|) shows that the field sample will have only 
^ 8% fewer group galaxies than a true global field sample. 
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We could correct for this bias by creating a true field sample 
which is an admixture made of 92% observed field sample 
and 8% group sample. However, all of our conclusions are 
insensitive to this correction, and for the sake of simplicity, 
we do not apply it to our results. 

Using these definitions yields a sample of 114 group 
galaxies and 128 field galaxies with 0.3 < z < 0.55 and My < 
-19. 



2.2 The z 0.1 Sample 

Our sample of low redshift galaxies is derived from the Mil- 
lennium Galaxy Catalog (MGC). The MGC is a 37.5 square 
degree B-band imaging survey carried out using the Wide 
Field Camera on the Isaac Newton Telescope (Li ske et aL 
I2OO3I ). The survey is a long, 35 arcmin wide strip, fully con- 
tained within both the Two Degree Field Galaxy Redshift 
Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS). 
The photometric catalogs are complete to B<24 mag, and 
the imaging is of sufficient quality to allow for the decom- 
position of galaxies into a bulge and disk component. 

MGCz, the redshift survey component of the MGC, was 
designed to obtain AAT/2dF spectra of B<20 gal axies which 
were not covered by either 2dFGRS or SDSS ( Driver et al.l 
I2OO5I ). This gives a redshift completeness of 96% for B<20 
galaxies. 



I2OO4I ). which is based on the 2dF survey, finds that ^ 55% 
of galaxies are in groups in the local universe. However, the 
fraction of 2PIGG groups with more than two members and 
within velocity dispersion limits of 100 km/s to 700 km/s 
(which are our criteria) is only 24%. The reduction of the 
sampling rate to match the CN0C2 completeness results in 
the non-detection of about half of these groups. A further 
30% of galaxies are not found in groups because the Carl- 
berg algorithm is not a strict friends-of-friends procedure: 
there is an additional step to check that candidate group 
galaxies are overdense with respect to the background. Ac- 
counting for these selection effects, we would only expect to 
find 9 % of galaxies satisfying our definition of a group. Fi- 
nally, we also remove groups which are not fully contained 
within the very narrow MGC strip, i.e. the group centers 
are within 500/if5^kpc of a survey edge. This step reduces 
the volume from which groups are selected from by ~ 30 
%. Our group catalogue is therefore incomplete relative to 
2PIGG, but our sample is robust (see ^4.3. ip and accurately 
reproduces the CN0C2 group-finding algorithm at higher 
redshift. The small number of group galaxies confirms, as 
suggested in ^2.1.1[ that our field sample is only slightly de- 
pleted in group galaxies when compared to the Universe as 
whole. Indeed, our results are unchanged if we include these 
group galaxies in our field population at this redshift, but 
for consistency with the CN0C2 groups, we do not. 



2.2.1 The Group and Field Samples 

There exist many low redshift group catalogues derived from 
either the SDSS or 2dF surveys. However, we wish to cre- 
ate a catalogue that can be compared directly and fairly 
with our higher-redshift group sample. The latter was de- 
rived from a two-step process — the initial survey and the 
targeted group follow-up — applied to an incomplete red- 
shift survey. This method is not the most direct or efficient 
way to find groups in our lower-redshift sample, but it does 
accurately reproduce our higher-redshift selection and the 
possible biases within. We discuss the method in detail in 
Appendix El Briefly, we reduce the completeness of the low 
redshift sample to match the sa mpling rate of the CN 0C2 
survey, and then reproduce the ICarlberg et ID (l200lh al - 
gorithm. To mimic our follow-up of Carlberg e t al.l ( 200ll ) 
groups, we then increase the completeness and recompute 
the group membership. 

Using this method we have found 19 groups with veloc- 
ity dispersions between 100 km/s and 700 km/s and which 
lie in the range 0.06 < z < 0.12. Using a volume- limited sam- 
ple of Mb < —18, we have a sample containing 99 group 
members and 3022 field galaxies. As discussed in § I2.1.1[ 
a true field sample would contain all group members and 
field galaxies. However, to maintain consistency between the 
CN0C2 and MGC samples we do not exclude the small num- 
ber of group galaxies from our field sample. In practice, due 
to the size of the field sample (~ 3000), adding or removing 
the 99 group galaxies has no effect on the bulk properties of 
the field. 

The low percentage of galaxies in groups in our sample 
seems in contradiction with other group catalogues based 
on local redshift surveys. For example, the 2d F Percola- 
tion Inferred Galaxy Groups (2PIGG) catalogue (|Eke et al.l 



2.3 Comparison of Surveys 

The measurement of a galaxy's morphology can depend on 
a number of factors besides its intrinsic morphology, such 
as the imaging wavelength, angular resolution and surface 
brightness limit. Because of these systematic differences, di- 
rect measurement of morphological evolution is difficult. In 
this paper, we largely concentrate on a direct comparison be- 
tween group and field galaxies at fixed redshift, which elimi- 
nates the effects of these systematic differences. Nonetheless, 
the MGC and CN0C2 are quite well matched in the key ar- 
eas, which allows us to compare the differential group and 
field behavior between redshifts. 

Specifically, the two surveys probe approximately the 
same restframe wavelength. The MGC survey is a rest frame 
B Band (observed frame ^ 440 nm) survey, while our ACS 
images are in the rest frame V Band (observed frame ^ 775 
nm). Since disks tend to be bluer than bulges one might 
expect a lower B/T when measured in the B-band; however 
the intrinsic morphological differences are only significant 
at much wider wavelength separation (|Tavlor-Mager et al.l 
I2OO7I ). 

The apparent surface brightness limit of the MGC sur- 
vey is 26 mag/arcsec^ and that of the HST ACS images is 30 
mag/arcsec^. Figure [2] shows the absolute surface brightness 
limits of the two surveys as a function of redshift, includ- 
ing (1+z)^ cosmological dimming. For the redshift ranges of 
interest, the absolute surface brightness limits are compara- 
ble. 

The excellent angular resolution of our HST ACS im- 
ages (point spread function (PSF) FWHM 0.1 arcsec) 
allows for morphological measurements of the CN0C2 sam- 
ple. In fact, this gives a physical resolution which is some- 
what better than the MGC survey (PSF FWHM 1 arcsec) 
in the redshift range of interest, as shown in Figure [3] Of 



Disks and Bulges of Group Galaxies 5 




2.5 



0.2 0.4 

Redshift 




0.2 0.4 

Redshift 



0.6 



Figure 2. A comparison of the absolute surface brightness hmits 
of our two samples of galaxies, as function of redshift. The solid 
line represents the MGC surface brightness (26 mag/arcsec^), 
while the dashed line is the equivalent for the CNOC2 sample 
(30 mag/arcsec^). 



Figure 3. A comparison of the point spread function FWHM of 
our two samples of galaxies as a function of redshift. The solid 
line represents the ground-based MGC PSF (1 arcsec), while the 
dashed line corresponds to the HST resolution of the higher red- 
shift CNOC2 sample (0.1 arcsec). 



greater concern is that the physical PSF of the MGC sur- 
vey has a large variation within the sample itself, due to 
the redshift range spanned by the galaxies. In Appendix [B1 
we investigate the effect of the physical resolution on the 
morphological measurements. We find that the resolution 
differences have no effect on the bulge-to-total light mea- 
surements, but do affect the asymmetry parameter, causing 
the measured asymmetry to be higher in galaxies with bet- 
ter resolution. In this paper, we only analyze asymmetry 
measures on matched samples of galaxies, which mitigates 
this effect. 

The propertie s of galaxies vary with luminosity (eg., 
iBaldrv et al.ll200 4) , so when comparing morphological prop- 
erties between the two surveys we must use a common lu- 
minosity limit. iFukugita et al have shown that late 
type galaxies (Scd to Sab) have B-V=0.5-0.78, while early 
type galaxies (SO and E) have B-V=0.85-0.96. Thus, our 
z = 0.4 magnitude limit of Mv<-19 corresponds approxi- 
mately to an MGC limit between Mb = —18.5 and —18.0, 
depending on type and neglecting any luminosity evolu- 
tion. Therefore, when directly comparing the group and field 
behavior between the two surveys we only consider MGC 
galaxies brighter than Mb ~ — 18 (except in Figure [5] which 
includes galaxies as faint as Mb = —16). 

2.4 Gim2d Morphological Measurements 

I Allen et"al] (|2006h have presented morphological measure- 
ments of the MGC and we follow their procedure to derive 
morphological parameters for the CN0 C2 sample. We use 
the parametric IRAF p ackage gim2d (|Simard et al.l I2OO2I : 
iMarleau Simar dl ll998l ). which fits the sky-subtracted sur- 
face brightness distribution of each galaxy with up to 
12 parameters describing a bulge and a disk component. 



GI m2d searches th e large -parameter space of models using 
a llMetropolis et aP (|l953l l algorithm, which i s inefficient but 
does not easily get trapped in local minima. iHaussler et al.l 
(2007) have shown that gim2d produces reliable fits with 
small systematic errors when the effective galaxy surface 
brightness is above the sky level, as it is for the galaxies in 
our sample. Q 

The gim2d algorithm can fit single component (Sersic 
profile) or two component (Sersic + exponential disk) mod- 
els to the galaxy profile. The Sersic profile is given by 

h{R) = /eexp [-6n[(i?/i?e)'/" - 1]] , (1) 

where le is the intensity at the radius. Re, and n is the 
Sersic index. The parameter bn is set to 1.9992n — 0.3271 
within gim2d to ensure that Re is the projected radius which 
encloses half the total luminosity. 

In a two component fit, the Sersic profile corresponds 
to the bulge model, and the disk is fit by an exponential 
model, 

Id{R)=IoeM-R/h), (2) 

where Iq is the central intensity, Id{R) is the disk light profile 
as a function of radius i?,and h is the scale length. 

Ideally, we would like to fit two components to all galax- 
ies, but often galaxies do not have two resolvable compo- 
nents. It is for this reason that we follow the prescription 

^ Haussler et al also show that an other parametric galaxy fitting 
code, GALFIT (|Peng et al.ll2002l ), performs better than GIm2d 
in crowded fields. However, GALFIT uses a downhill gradient 
algorithm which, although very efficient, may not be as robust 
as the simulated annealing technique of GIm2d for problems with 
many local minima, like two- component fitting. 
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of [All en et alJ (|2006l ). who have made a careful study of 
GIm2d output. They suggest a logical filtering system which 
initially fits a two-component model to the galaxy. For galax- 
ies which have "normal" light profiles, this fit is kept, but 
for those galaxies which have perturbed profiles or are obvi- 
ously better described by a single component, we fit a pure 
Sersic profile. In practice, this has little effect on our results. 



2.4-1 Residual Substructure 

To quantify the substructure in the surface brightness 
profiles we u s e the residual parameter, R, as defined by 
ISchade et all (|l995l l. This parameter is also known as the 
"asymmetry parameter, R" (|lm et a l. 2002) and the "resid- 
ual substructure parameter, S" (jMcIntosh et al.ll2004 V It is 
defined as 



R = Rt + Ra 



with 



Rt 



Ra 



S(|i?.,+i?Jf 1/2) 



1/2) 



(3) 



(4) 



(5) 



where Rt is the total residual parameter and Ra is the 
asymmetric residual parameter. Kij is the flux at pixel po- 
sition (i,j) in the residual image, ^jj^ is the flux at (i,j) in 
the residual image rotated by ISO''. Bij and Bjf^ are the 
corresponding values in the background noise. Finally, lij is 
defined as the flux at (i,j) in the object image. The sum is 
done over all pixels out to r= 2r hi, where ihi is the radius 
at which half the galaxy's light is enclosed. 



3 RESULTS 

Although our data sample is comprised of a relatively large 
number of group galaxies (99 in the MGC sample and 
114 in the CN0C2 sample), each group has typically less 
than ten spectroscopically conflrmed members. Therefore, 
we stack the individual groups to maximize the s ignal of the 
bulk group galaxy proprieties. IWeinmann et al.l (|2006l l have 
shown that the velocity dispersion of a group is a poor tracer 
of the mass of the group halo, especially for groups contain- 
ing few confirmed members. Although our group velocity 
dispe rsions vary from 100 km/s to 700 km/s. 1wilman et al.l 
(|200 5b) has shown that the individual CN0C2 groups show 
no significant differences based on group type or velocity 
dispersion from a combined group. Thus we combine all the 
galaxies within SOO/ifg^kpc of a group center at each redshift 
to form a stacked Z'-^O.l group and a stacked z^O.4 group. 
Because of the large uncertainties on the velocity disper- 
sions we do not attempt to estimate a "virial radius" for 
each group, but instead simply require a group member to 
be within SOO/ifg^kpc of the group center, corresponding ap- 
proximately to the expected virial radius of a typical group 
in our sample, with a 360km/s velocity dispersion. 



3.1 B/T Distribution 

In Figure S] we present the quantitative morphology dis- 
tribution of the two samples. We compare field and group 
distributions of the ratio of the bulge luminosity to total 
luminosity (B/T) of each galaxy. Pure bulge galaxies have 
a B/T of 1, while pure disk galaxies have a B/T of 0. The 
left hand panel of Figure H] shows the B/T distribution for 
the z ^ 0.4 sample of galaxies. The red line represents the 
stacked group from the CN0C2 sample and the black line 
represents the field galaxies. A Kolmogorov-Smirnov (KS) 
test shows that the group and field sample are not drawn 
from a common parent distribution, with 98.4% confidence. 
The most significant difference within the CN0C2 sample 
is at B/T<0.2, where the fraction of galaxies in the field 
(-- 57 ± 5%) is higher than that in groups 41 zb 5%). 
In the right hand panel of Figure H] we present the B/T 
distribution for the z ^ 0.1 sample. Similar to the CN0C2 
sample, the fraction of B/T < 0.2 galaxies is much higher 
in the field (- 54 zb 1%) than in the groups (- 32 ± 6%). A 
KS test rules out a common origin for the group and field 
distributions of the MGC sample at greater than 99.9 % 
confidence. 

In both plots we have seen evidence for the well-known 
morphology-density relation. In this case, it is manifested 
as a deficit of disk dominated galaxies in the group sam- 
ples when compared to the field. However, it is known that 
bright galaxies tend to be more frequently bulge-dominated 
than faint galaxies. Therefore, it is possible that this form of 
the morphology-density relation is related to different field 
and group luminosity distributions rather than any intrinsic 
difference between group and field galaxies at fixed lumi- 
nosity. To explore this we present Figure O which shows 
the fraction of disk-dominated galaxies as a function of to- 
tal galaxy magnitude. Studies of the B/T distribution of 
visually classified galaxies have shown that elliptical and SO 
gala xies predominately have B/T>0.3-0.4 (i Tran et al.]l200ll 
and lWilman et al.ll2008l ). but we chose to define "disk dom- 
inated" to indicate galaxies with B/T<0.2. This choice is 
made to isolate those galaxies in the bin with the most sig- 
nificant difference between the group and field samples, but 
our conclusions are unchanged even if we use B/T<0.4 to 
define disk-dominated galaxies. 

The left hand panel of Figure [5] shows the z ^ 0.4 sam- 
ple, with the thick black line representing field galaxies and 
the thin red line representing the groups. In any one luminos- 
ity bin there is no significant difference (ie.> 1 a) between 
the fraction of disk galaxies in the group and field; how- 
ever, overall there is a systematic difference of 5.5 ± 2 %, 
in the sense that the fraction of disk-dominated galaxies in 
groups is always lower than in field galaxies of comparable 
luminosity. 

The right hand panel of Figure [5] shows the disk fraction 
as a function of magnitude for the z ^ 0.1 sample. The 
magnitude is the rest frame B band, the waveband in which 
the B/T decomposition was done. In this figure, we show 
all galaxies in our redshift range to Mb=-16, fainter than 
our volume- limited sample of Mb<-IS. The dashed vertical 
line indicates the equivalent B-band limit of the CN0C2 
Mv = — 19 magnitude limit, assumin g a B-V color of 0.63 , 
which is typical of a late- type galaxy ([Fukugita et al.|[l995l ) . 
At this redshift there is clearly a significant difference in the 
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Figure 4. The histogram of relative bulge luminosity for the z ~ 0.4 (left) and z ~ 0.1 (right) samples. The thin red line in each plot 
indicates the stacked group at that redshift, and the thick black line is the field sample. Uncertainties are estimated with the jackknife 
technique. 




Figure 5. The fraction of galaxies with B/T<0.2 (disk-dominated) as a function of absolute magnitude. Left: The z ~ 0.4 sample. 
Magnitudes are measured using the flux calculated by GIm2d and k-corrected to rest frame V band. Right: The z ~ 0.1 sample, where 
magnitudes are measured in the B band. The fleld galaxies are represented by the thick black line and the group galaxies by the thin red 
line. The dashed vertical line indicates the equivalent B-band limit corresponding to the CNOC2 My = — 19 magnitude limit, assuming 
a B-V color of 0.63. Uncertainties are measured using the jackknife technique. 



fraction of disk- dominated systems between group and field 
galaxies of the same luminosity. This difference is 24 it 6 % 
over the full magnitude range of the sample, and 19 ± 6 % 
brighter than the equivalent CN0C2 luminosity limit. 

Although the B/T distributions of the group and field 
galaxies are different at both redshifts, this predominately 
reflects a difference in luminosity distributions at high red- 
shifts, but an intrinsic difference in the fraction of disk galax- 
ies of fixed luminosity at low redshift. It therefore appears 



that the morphological segregation in groups has increased 
significantly from z~0.4toz'-^0.1. We recall that the z=0.1 
sample is measured in the rest frame B band while the z=0.4 
sample is measured in the rest frame V band. It would be 
possible to mimic our results if the disks of group galax- 
ies were significantly redder than the disks of field galaxies. 
However, to create a B-band disk fraction consistent with 
the V band disk fraction, the disks of group galaxies must 
be B-V=2.07 redder than the disks of field galaxies. On aver- 
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age, in the MGC sample, pure disk (B/T=0) group galaxies 
are only B-V=0.01 redder than the pure disk field galaxies. 
Thus, it appears the evolution is real. We will further ex- 
plore the colours of these galaxies, including a comparison 
with the CN0C2 sample, in a future paper. 

3.2 Structural Parameters 

The evolution in morphological segregation suggests there 
is a change occurring within the group galaxies between z 
'^0.4 and z ~ 0.1. Here we focus on the possible causes 
of this evolution. We first look for structural differences be- 
tween the group and field galaxies. If star formation in group 
galaxies is quenched, then the disk of the galaxy may slowly 
fade away, and we might expect to see a departure from the 
normal scaling relations between disk size and luminosity. 

In Figure [6] we plot the distribution of disk magnitudes 
as a function of the disk scalelength, excluding only the most 
bulge-dominated galaxies (B/T > O.7).0The left hand panel 
shows the z ^ OA sample, while the right hand panel shows 
the z ^ 0.1 sample. In both plots it is evident that there 
is a correlation between the disk size and its brightness, 
such that brighter disks have larger disk scale lengths. This 
is not surprising, but the fact that the group and field lie 
on the same relation (although with large scatter) is. The 
thick solid black line is the best fit to the scaling relation 
of the field galaxies. Best fit lines are determined using a 
robust biweight estimator which minimizes the effect of dis- 
tant outliers. In principle, uniform fading of a perfectly ex- 
ponential disk would result in a lower disk luminosity, but 
would leave the scalelength unchanged Therefore, an ideal 
population of faded disk galaxies would exhibit the same 
scaling relation, but with different normalization. Adopting 
this assumption, we fit only the normalization to the group 
galaxy relation, while maintaining the slope defined by the 
field galaxies, as shown by the thin red line. There is no sta- 
tistically significant difference between the normalizations 
of the group and field populations (a difference of 0.086 =b 
0.134 for the z=0.4 sample, and 0.103 ± 0.151 for the z=0.1 
sample) . 

Since we do not see any significant difference in the disk 
scaling relations for group galaxies, it may be that the pro- 
cess of morphological transformation is instead dominated 
by a growing bulge (for example, through mergers). Again, 
if this were true, we might expect to see a deviation in the 
group and field bulge scaling relations. In Figure [71 we plot 
the distribution of bulge magnitudes as a function of the 
bulge half light radius, excluding the most disk-dominated 
galaxies (B/T < 0.3). The left hand panel shows the z ^ 
0.4 sample and the right hand panel has the z^O.l sample. 

^ An occasional problem in modeling the surface brightness pro- 
file of galaxies with automated programs is the tendency to fit 
small disk (or bulge) components to galaxies which may have 
(e.g.) twisted isophotes. Thus we restrict our analysis to galaxies 
which have significant disk (or bulge) components when looking 
a t their scaling p roperties. 

^ iHaussler et al.l (|2007l ) have shown that the recovered scale 
length is underestimated for low S/N galaxies which are fit with 
a Sersic profile. However, our data are sufficiently deep to avoid 
this problem, typically reaching the surface brightness limiting 
isophotal radius at 3-5 disk scale lengths. 



The field and group distributions are again similar, at both 
redshifts. In 3H we examine the constraints these findings 
place on the amount of fading which is possible in the group 
sample. 

3.3 Asymmetry 

In this section we examine the asymmetries of the galaxies, 
to try to untangle a merger-driven transformation scenario 
from a gradual disk fading model. Galaxy mergers and ha- 
rassment often produce noticeable asymmetric features, like 
tidal tails, which could manifest themselves as deviations 
from the smooth surface brightness profiles we have used 
for the morphological measurements. If galaxy mergers were 
enhanced in groups, we would expect to see an increase in 
the fraction of galaxies with large asymmetries. On the other 
hand, the cessation of star formation in the disk would likely 
result in the disappearance of bright clumps and spiral arms, 
thereby removing residual substructure, making the disks 
appear smoother. 

We have seen that the fraction of disk galaxies depends 
on magnitude, and that the difference in the high redshift 
B/T distribution is partly due to different luminosity distri- 
butions in the field and group samples. For this reason, in the 
rest of this section, we match the field and group luminos- 
ity and redshift distributions. In the MGC sample, we match 
each group galaxy to the field galaxy with the closest magni- 
tude and within 0.03 in redshift. The magnitude differences 
for these matched galaxies are all less than 0.02 because of 
the large number of available field galaxies in the sample. 
Similarly, for the CN0C2 sample, we match group galaxies 
to the nearest field galaxy in magnitude (within 0.2 mags) 
and within 0.05 in redshift. We match all 99 group galaxies 
in the MGC and 105 of the 114 CN0C2 group galaxies. The 
brightest CN0C2 group galaxies have no field counterpart, 
and so are effectively excluded from the remainder of this 
analysis. 

To probe the substructure of the galaxies, we have cal- 
culated the asymmetry parameter according to the defini- 
tions given in ^2.4.11 In Figure [8] we show the distribution in 
asymmetry for bulge dominated (B/T>0.5) and disk dom- 
inated (B/T<0.5) galaxies for both the CN0C2 and MGC 
surveys. The dashed vertical line indicates the lower limit for 
"highly asymmetric" galaxies {Rt + Ra > 0.16), as defined 
by Pat ton et al. (2005 ). Disk dominated galaxies have much 
higher median asymmetries than bulge dominated galaxies, 
as expected. However, there is no appreciable difference in 
the median asymmetry between the group and luminosity- 
matched field, in any of the samples. Although we can not 
resolve the disapperance of individual HII regions, visual in- 
spection of disk dominated galaxies in Figures ICll and IC2I 
clearly show that galaxies with low B/T and high asym- 
metries have strong asymmetric structures typical of star 
forming galaxies, ie. large spiral arms and lumpy regions. 
The lack of a systematic difference in the asymmetries of 
matched group and field disk-dominated galaxies means that 
we find no evidence for a mechanism that suppresses star 
formation in group disks. 

We note that, as shown in Figure O the CN0C2 sample 
has a lower physical PSF than the MGC sample. We show in 
Appendix [B] that galaxies which are blurred to have a lower 
physical resolution have lower asymmetry values. Therefore, 
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Figure 6. The distribution of disk magnitudes as a function of disk scalelength for non-bulge dominated galaxies (B/T <0.7). Left: z ~ 
0.4 sample. Right: z~ 0.1 sample. Large red points are group galaxies and small, black points are field galaxies. The thick solid black line 
is the best fit to the field, thin red line is the group best fit. The dotted (dashed) black line shows the effect of truncating star formation 
in field galaxies 1(3) Gyrs ago. 
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Figure 7. The distribution of bulge magnitude as a function of bulge half light radius for non-disk dominated (B/T > 0.3) galaxies. 
Left: The z~ 0.4 sample. Large, red points are group galaxies and small, black points are field galaxies. Right: The z~ 0.1 sample. Red, 
filled points are group galaxies and black points are field galaxies. 



we do not consider the change in median asymmetry between 
z ^ OA and z ^ 0.1 as evidence of real evolution. 

4 DISCUSSION 

4.1 Galaxies in Transformation? 

In this section we investigate how our results in 33] might 
constrain the number galaxies that are in the process of 



transforming within the group environment. Specifically, we 
have shown in ^3.21 that the bulge- and disk-components 
appear to obey scaling relations between size and luminos- 
ity that are independent of environment. This suggests that 
any transformation mechanism must either leave these scal- 
ing laws intact, or only affect a small number of galaxies at 
the epoc h of observation. T o inves tigate this, we consider 
a simple iBruzual k CharlotI (|2003h model with a constant 
star formation rate as appropriate for the disk components. 
If this star formation is suddenly truncated, within 1 Gyr it 
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Figure 8. The distribution of asymmetries (Rt+Ra) for CNOC2 
bulge dominated galaxies (top left), CNOC2 disk dominated 
galaxies (top right), MGC bulge dominated galaxies (bottom 
left), and MGC disk dominated galaxies (bottom right). In all 
plots group galaxies are represented by the thin red line and the 
luminosity matched field is the thick black line. The dashed verti- 
cal line indicates the lower limit fo r "highly asymmetri c" galaxies 
{Rt ^ Ra> 0.16) as defined bv lPatton et al.l (|2005[ ). 



(4 zb 2 %) for the MGC sample. The upper limit for the 
z—OA sample is quite high, as expected given the relatively 
small sample size. Therefore, based on these data alone, we 
cannot rule out the hypothesis that a substantial fraction of 
these galaxies are undergoing a significant change in their 
star formation rate. The upper limit for the z = 0.1 sam- 
ple is much more restrictive, partly because of the sample 
size and partly because the rest-frame B-band is more sen- 
sitive to recent star formation. Our limits mean that any 
mechanism able to truncate star formation in disks is not 
dominant in present day groups. 

Mergers are likely to be rare, but transformative events, 
so an increased merger history in groups would not be ex- 
pected to increase the median asymmetry of group galaxies, 
but rather increase the fraction of highly asymmetric galax- 
ies. By examining the fraction of galaxies with Rt + Ra > 
0.16, we see that an extra 6 =b 3 % of group galaxies are 
highly asymmetric in the bulge-dominated samples at both 
redshifts, compared with the matched field sample {|8]). The 
fraction of highly asymmetric galaxies in the disk-dominated 
samples is similar for both the group and the luminosity- 
matched field. This may indicate that there is an increase 
of merging or interacting galaxies in groups. A visual in- 
spection of the CN0C2 galaxies shows that 4 out of 9 of 
the "high- asymmetry" bulge dominated galaxies are indeed 
merging or interacting, as shown in Appendix [Cl 



will have subsequently faded by 1.23 magnitudes in lest-B 
(as measured for the MGC sample), or 0.68 magnitudes in 
lesi-V (appropriate for the higher-z CN0C2 sample). Af- 
ter 3 Gyr, the amount of fading expected is 1.91 mags (B) 
or 1.39 mags (V). However, such truncation should not af- 
fect the measured scale length of the disk. Therefore, if star 
formation were truncated in the entire population of group 
galaxies, we would expect the normalization of the scaling 
relation to change by these amounts. These relations are 
shown by the dotted (1 Gyr) and dashed (3 Gyr) black lines 
in Figure [6] These lines are significantly offset from the mea- 
sured group relation (which is consistent with that of the 
field); therefore, we can easily rule out that star formation 
has been recently truncated in the entire disk population. 

We would next like to establish what fraction of group 
galaxies could have undergone 1 (3) Gyr fading and still 
have a scaling relation that is consistent with the observed 
field scaling relation (ie. < 2 a difference in the normaliza- 
tion, corresponding to a difference of 0.2 magnitudes for the 
MGC sample, and 0.18 magnitudes for CN0C2). To assess 
this we randomly choose a sample of group galaxies from 
Figure [6] and fade their disks for 1 (3) Gyrs according to the 
Bruzual & Chariot model described above. We recompute 
the B/T ratios of the sample, and impose our selection cri- 
teria on total magnitude and B/T (recall we are excluding 
the most bulge-dominated galaxies from this analysis). The 
normalization of the scaling relation is then refit to the new 
set of data, keeping the slope fixed. We can then determine 
an upper limit on the fraction of galaxies which may have 
undergone such fading for 1 (3) Gyrs; these limits are 41 
± 3 % (29 ± 4 %) for the CN0C2 sample, and 9 ± 3 % 



4.2 A comparison with X-ray selected groups 

In ^3.3[ we have shown that the median asymmetry of 
group and field galaxies are statisti cally indist inguishable. 
These results are interesting because Tran et al.l (j2001) have 
shown that, in a sample of local X-ray selected groups, there 
is evidence for smoother disks in group galaxies than in 
the field. Studies of blue cluster galaxies have also shown 
that they have significantly lo wer asymmetry value s than 
their blue field counterparts ("Mcintosh et al.^ |2004| V Per- 
haps most intriguingly, Homeier et al. (2006) has recently 
shown that X-ray luminous clusters have galaxies with sig- 
nificantly lower average asymmetries than X-ray faint clus- 
ters. A related difference from our results comes from stud- 
ies of X-ray group galaxy morphology, which have sh own 
that the fraction of early types is ^ 0.7 ([^^Ichaev et al.l 
I2OO6I : IJeltema et al.1 [20071 ) in the same magnitude and red- 
shift range as our CN0C2 sample. To compare our data 
with this nu mber we define ear ly-type galaxies as those with 
B/T> 0.4, as lTran etld] (|200l[ ) has shown that this provides 
a good match with early type galaxies as classified on the 
Hubble-sequence. Using this definition, we find an early-type 
fraction of only 46 % in the CN0C2 group sample, con- 
siderably lower than found in X-ray groups at this redshift. 
We will revisit this issue in a future paper when we consider 
and compare the results of Hubb le-sequence morphological 
classification (|Wilman et ahllioosh . These results, combined 
with our result that the median asymmetry in optically se- 
lected group galaxies is not different from the luminosity 
matched field, point to the role that the hot intergalactic 
medium (IGM) may play in the smoothing of disk galaxies. 
Alternatively, the progenitors of optically selected groups 
may be different from the progenitors of X-ray selected ones. 
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Figure 9. The distribution of halo masses for each galaxy in the I Bower et al.l (|2006|) model catalogue (thin, black line) and the 
distributions of galaxy halos in the simulated groups (thick, red line). The high redshift sample is on the left and the low redshift sample 
is on the right. 



4.3 Comparison with semi-analytic galaxy models 

Recently, large dark matter simulations of large volumes 
have allowed theorists to produce usefully large catalogues 
of model galaxies, employing detailed modeling of galaxy 
formation based on relatively simple prescriptions for rel- 
evant physical processes. In this section we use the cata- 
logues of one such "semi-analytic" galaxy formation model 
(|Bower et al ] l2Q06 h to compare with our data. The Bower 
et al. model uses th e dark matter Millennium simulation 
(jSpringel et al.ll2005l ). a ACDM cosmological box with 500//i 
Mpc sides, as the basis for the merger trees. The algorithm 
is based o n the earlier GALFORM models of Benso n~et al.l 
(2003) and Cole etHI (jioO O). The principal change from the 
Benson et al. model is a prescription for the quenching of star 
formation in massive halos by feedback powered by accre- 
tion onto a supermassive black hole. However, perhaps the 
most important change for our purposes is the modification 
of the method for computing disk instabilities. Disk instabil- 
ities are now the dominant mode of bulge formation in these 
models, although the brightest galaxies are still more often 
formed through mergers. Unlike older models, morphology 
is now sensitive to the baryonic physics of disks, rather than 
the (more robustly-predicted) merger history. 



4.3.1 Constructing a mock catalogue 

We construct mock catalogues using the lBower et al.1 (|2006l l 
model for the z=0.1 and z—OA redshift time steps. Because 
of the large size of the Millennium simulation and the small 
spread in redshift of each of our samples, a single simula- 
tion box can be used for each epoch. Using the same group- 
finding procedure as described in detail in Appendix [A] we 
construct a stacked group sample to compare directly with 
our observations. 

In Figure [9l we present the resulting distribution of 



galaxy halos in our group sample, compared with all halos in 
the Millennium simulation. We see that at both redshifts the 
group-finding algorithm selects predominately galaxies that 
are in halos with masses 5 x 10^^ < Mhaio/M© < 10^^. Both 
samples peak at Mhaio ^ 6 x lO^^M©, but the CN0C2 sam- 
ple is somewhat biased toward lower masses than the MGC 
sample. Our algorithm finds a large percentage of groups 
that are made up principally of large dark matter halos: 
76.4 % (85.1 %) of z=0.4 (z=0.1) galaxies are in halos with 
Mhaio > 10^^ M0. In both plots there is a second peak - a 
distribution of low mass halos - which are contamination. 
However, we find that 79% (z=0.4) and 89% (z=0) of the 
Mhaio <10^^ M0 galaxies are within bOOh^^kpc of a galaxy 
which resides in a halo with Mhaio >10^^ M©. Thus, the ma- 
jority of our "contamination" is due to galaxies on the out- 
skirts of a true group. This confirms that the CarlbeT g et al.l 
algorithm selects groups which are real, and repre- 
sentative of massive dark matter halos, as also confirmed 
by previous weak lensing measu rements (|Parker et al] 
and our follow-up spectroscopy (Wilman et al. 2005b). Only 
2.5 % (3.1 %) of our z=0.4 (z=0.1) galaxies are not associ- 
ated with a massive (Mhaio > 10^^ M©) halo. 

4.3.2 B/T distribution 

In Figure [TOl we show the 'Bower et al.l (|2006 h model predic- 
tions for the B/T distributions corresponding to our data 
samples. In all four panels the data are shown with a thick 
black line and the model predictions are shown with a thin, 
dashed, red line. The models are limited at Mb < — 18 in 
the MGC comparisons, and by < -19 for the CN0C2 
comparison. 

It is clear that there is remarkable agreement between 
the model and the data at z=0.1, especially considering that 
the model does not take into account any observational un- 
certainties associated with deriving B/T ratios from the sur- 
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Figure 10. The B/T distributions of the iBower et all (|2006l ) 

model, compared with the observed group and field samples at 
z=0.1 and z=OA. The model is the thin, red line in each panel 
and the data is the thick black line. 
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Figure 11. The relative disk deficiency in groups, parametrized 
as 1-group disk fraction/field disk fraction, where a disk has a 
B/T < 0.2. This is shown at z=0.4 (right panel) and z=0.1 (left 
panel) for both the Bower et al model (thin red line) and the data 
(thick black line). 



face brightness profile alone. However, the agreement be- 
tween the models and the data at z=0.4 is not as good. 
In particular, the model underpredicts the fraction of B/T 
< 0.2 galaxies and overpredicts the fraction of B/T > 0.8 
galaxies in both the groups and the field. Intriguingly, Fig- 
ure [10] shows the models also predict that the fraction of 
disk dominated galaxies increases in the field between z=0.4 
and z=0.1, but remains constant within the groups. 

Given these predictions, we are now encouraged to in- 
vestigate the time evolution of the disk fraction as a function 
of luminosity. From our data we have seen that the differ- 
ences in the disk fraction as a function of magnitude are 
small at z=0.4 5.5 di 2 %), but quite large in the local 
universe (~ 19 it 6 %). To address this, we present Figure 
llli which shows the relative disk (B/T < 0.2) deficiency 
between the field and the group samples, as a function of 
luminosity. The disk deficiency is the difference between the 
group and field disk fraction divided by the field disk frac- 
tion. This gives a measure of the fraction of field disks which 
are absent at a similar magnitude in the groups. In the left 
panel of Figure [TT] we show the low redshift sample. The 
data from the MGC sample agrees well on average with the 
model predictions. In the right panel of Figure [TT] we show 
the same comparison but for the high redshift sample. Al- 
though the average value of the group disk deficit is correctly 
predicted at z=0.1, the Bower et al. model predicts a group 
disk deficit which is much higher than the data at z=0.4. We 
note that, in the model, the predicted evolution is of a sim- 
ilar magnitude whether measured consistently in rest frame 
V or rest frame B, indirectly supporting our argument that 
the observed evolution is not driven by the difference in rest 
wavelength sampled by the two surveys. 

We have seen in Figure [TOlthat the models underpredict 
the fraction of B/T < 0.2 at z=0.4 and that this leads to an 
overprediction of the disk deficiency in groups at the high 



redshift epoch. Intriguingly, inspection of Figure \T0\ shows 
that this is because the models predict that the fraction 
of disk dominated galaxies increases in the field between 
z=0.4 and z=0.1, but remains constant within the groups. 
While a direct comparison between the two epochs could be 
complicated by the different observed wavebands, we note 
that the same model predictions exist when considering only 
the underlying stellar mass. This is at odds with the data, 
which show that the B/T < 0.2 fraction decreases in the 
groups, and remains constant in the field. 



5 CONCLUSIONS 

We have presented a quantitative morphology study of 
optically-selected galaxy groups from two redshift surveys: 
CN0C2 (z'-^0.4), supplemented with significant additional 
Magellan spectroscopy, and MGC (z^O.l). We have com- 
pared these data with a similarly selected sample of groups 
drawn from the se mi-analytic galaxy formation models of 
iBower et al] (|2006l ). Our findings are: 

• There is a significant difference, as indicated by a KS 
test, in the fractional bulge luminosity (B/T) distribution 
of group and field galaxies, in both the high and low red- 
shift samples. The dominant difference is the deficit of disk- 
dominated (B/T < 0.2) galaxies in the group samples. 

• The difference in the disk fraction (B/T<0.2) of group 
galaxies relative to the field shows significant evolution be- 
tween z=0.4 and z=0.1. At a given luminosity in the CN0C2 
sample, the groups have ^ 5.5 ± 2 % fewer B/T<0.2 galaxies 
than the field. By z=0.1 this difference has increased signif- 
icantly, so that groups have ~ 19 =b 6 % fewer B/T<0.2 
galaxies than the field in the same magnitude range. Al- 
though the z=0.1 sample traces rest frame B while the z—OA 
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sample traces rest frame V, this difference is unlikely to be 
able to explain the differences. 

• At neither redshift do we see any evidence that the bulk 
properties of the existing disks are significantly different for 
group galaxies than for field galaxies. They lie on similar 
scaling relations and show similar asymmetry distributions. 
There is no evidence that groups are actively perturbing or 
otherwise affecting a large fraction of the group disk popu- 
lation. 

• We find that there is a small enhancement in the 
fraction of bulge- dominated group galaxies that are highly 
asymmetric, relative to bulge-dominated galaxies in the 
field. This may be consistent with enhanced merging in the 
group environment. Visual inspection of high asymmetry, 
bulge-dominated CN0C2 galaxies shows that 44% (4/9) ex- 
hibit clear evidence of interactions or merging. 

• A sample of galaxies drawn from the semi-analytic 
galaxy catalogues of iBower et al.l (|2006l ) was shown to agree 
remarkably well with the B/T distribution of the field and 
group galaxies at z=0.1. However, our data have shown that 
time evolution of the B/T distributions predicted by the 
models is not seen in our data. In particular, the Bower 
et al. model underpredicts the fraction of disk dominated 
galaxies at z ^ 0.4. 

The morphological difference between group and field 
galaxies at z=0.4 is mostly due to the tendency for group 
galaxies to be more luminous and, therefore, more bulge- 
dominated than field galaxies. This is consistent with our 
previous findings about group galaxies; namely, that their 
M/L ratios are consistent w ith the passive evolu tion of a pre- 
dominately old population tealog h et ahlbOOTl l , because the 
dominant difference in group galaxies is their pre-existing 
tendency to be bulge-dominated. This is also consistent with 
the fairly small difference in the emission-line fraction of 
group and field gala xies, at fixed stellar rn ass (|Balogh et al.l 
I2OO7I ) or magnitude (jWilman et"alll2005al). 

The failure of the iBower et al.l (|2006l l to reproduce 
the time evolution of the group disk deficit is interest- 
ing. These type of models predict a fairly rapid "stran- 
gulation" of galaxies once they enter larger halos, which 
causes the star formation rate to decrease on timescales of 
a few Gyr. This starts to play a large role at ~10^^ solar 
masses, significantly less massive than most of our groups. 
.Weinm ann et al. (2 006 ) have shown that this mechanism 
is too effective, and produces a homogeneous, red satellite 
population at all magnit udes, in groups and clusters, which 
is not observed at z=0. iGilb ank Baloghl i 20081 ) have re- 
cently shown that this problem extends out to at least z=l. 
It is likely that the incorrect disk deficit evolution is an- 
other manifestation of the maximally efficient star formation 
quenching mechanism used by the models. 

We have found that, in contrast to X-ray selected 
groups, our optically-selected group galaxies have median 
asymmetries that are similar to field galaxies. This may 
point to a possible role for galaxy interactions with the 
hot IGM. Alternatively, the progenitors of X-ray selected 
group galaxies may be fundamentally different from the 
progenitors of optically selected group galaxies. X-ray se- 
lected groups are more likely to be relaxed, virialized struc- 
tures, suggesting that they were assembled earlier than op- 
tically selected groups. In the current models, the stran- 



gulation mechanism is assumed to operate efficiently in 
small haloes. Thus, little difference is expected in the pop- 
ulation of different types of groups many Gyr later (i.e. 
at the epochs of interest here), since star formation has 
long ceased in most group members. However, recently it 
has been suggested that infalling galaxies may be able to 
retain a significant fraction of t heir gas (|McCarthv et al.] 
I2OO7I : 'Kawata &: MulchaevI l2007l l, significantly increasing 
the timescale for star formation to decrease. In this case, 
the headstart given to galaxies that fall into groups a little 
earlier may be better able to explain the difference between 
X-ray and optically selected groups. 
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APPENDIX A: LOW REDSHIFT GROUP 
FINDING ALGORITHM 

The primary goal of our low redshift group finding algorithm 
is to reproduce the selection criteria applied to our high red- 
shift groups. Thus, our algorithm is not the most efficient 
method possible, but it does accurately reproduce our se- 
lection and the possible biases within. Further, this method 
will not result in a complete sample of groups in the MGC 
strip. 

We first find groups in the SPSS main galaxy sample us- 
ing the original method of ICarlberg et al.1 (|200lh . The MGC 
strip is a narrow region, ~ 35 arcmins across, which makes 
group finding within the strip itself difficult. For this reason, 
we first find groups in the SDSS, in a region 2 degrees across 
and centered on the MGC. We define our SDSS galaxy sam- 
ple to be directly analogous to the CN0C2 sample. There 
are two areas of particular relevance to this work where 
these differ: completeness and depth. Because the SDSS has 
a much higher completeness 90%) than the CN0C2 red- 
shift survey (~ 48 %), we randomly remove half the SDSS 
galaxies. Further, we use the same absolute magnitude cut 
as Carlberg et al., Mr = —18.5, with an additional evolution 
correction of 1 magnitude per unit redshift. 

Carlberg et al.'s primary goal was to find virialized 
groups in overdense environments, so they estimate the over- 
density of each galaxy and restrict their group finding al- 
gorithm to galaxies in dense environments. A cylinder of 
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Figure Al. The redshift distribution of the full sample of SDSS 
galaxies which were used to find groups at low redshift. The red 
line is the Maxwellian fit to the data which was used to estimate 
the background. 

0.33/1^5^ Mpc radius and ±6.67h^^ Mpc line-of-sight depth 
is centered around each galaxy and the number of galax- 
ies within the cylinder is counted. If there are fewer than 
3 neighbors in this cylinder, the process is repeated with 
a cylinder 1.5 times larger. A background estimate is then 
obtained by randomly drawing points from a redshift distri- 
bution fit to the entire sample. Figure [All shows the redshift 
distribution of galaxies in our sample and our analytic fit. 
If the number of neighbours in the cylinder is greater than 
the background estimate then the main galaxy is kept as a 
possible group member. 

Starting with the galaxy with the greatest overdensity, 
we begin a trial group by adding any galaxies within the 
original cylinder, and any of their friends. When we run out 
of friends we have a trial group, for which the geometric posi- 
tion, redshift and velocity dispersion are computed. Galaxies 
are trimmed or added within 1.5i?2oo and three velocity dis- 
persions, where R200 is the radius at which the density is 200 
times the critical density. This process is iterated four times 
with the requirement that the last two iterations are identi- 
cal. A group is moved to the next stage if it has more than 
two members. This concludes the Carlberg et al. algorithm. 

The nex t stage is to emulate the process in 
I Wilman et al.l (|2005b[ ), to account for the targeted spectro- 
scopic follow-up, which resulted in a more nearly complete 
redshift sampling around selected groups. We do this by in- 
cluding the complete SDSS catalogue. We use the Carlberg 
initial centres but set the velocity dispersion equal to 500 
km/s, as was done by Wilman et al. This was done to remove 
any bias in the starting velocity dispersions, which were only 
based on very few galaxies. We again iterate on these po- 
sitions using the entire SDSS catalogue and recompute the 
luminosity weighted centres. We compute the velocity dis- 
persion at each step using the Capper estimator and remove 
galaxies outside two velocity dispersions and 500/1^5^ kpc. Fi- 
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nally, we keep only those groups which he completely within 
the MGC strip. Using this method we have 19 groups with 
velocity dispersions between 100 km/s and 700 km/s, and 
which lie within 0.04 < z < 0.12. 

There exist a large number of group and cluster cata- 
logues based on the SDSS and 2DF surveys with which we 
can compare. This is especially important to calibrate the 
systematic effects which may be present in our high red- 
shift sample, which doesn't have sufficient completeness to 
quantify within the survey itself. One of the more popu- 
lar group finding algorithms, and the rnost di rect analogue 
to our method, is that of Berl ind et all (|2006l l. They use a 
traditional friends-of-friends algorithm in position-redshift 
space to find groups in three different volume limited sam- 
ples. We find that 12 out of 13 of our groups below z=0.1, 
the depth of the deepest Berlind sample, are also found in 
the BerUnd catalogue, ie. the Berlind group centres are con- 
tained within our groups. Eighty-five of the 99 group mem- 
bers which make up the low redshift group sample in this 
paper would also be group members if we were to use the 
Berlind catalogue as our group catalogue. 



APPENDIX B: DEPENDENCE OF GALAXY 
PROPERTIES ON PSF SIZE 
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As discussed in ^TSl the CN0C2 and MGC surveys com- 
pare well in absolute surface brightness limits, physical size 
of the PSF and the rest waveband used in the morphological 
decomposition. In fact, the biggest variation in these param- 
eters is actually within the MGC sample itself. The PSF size 
is ~ 3 kpc in size at z = 0.14 and just ~ 1 kpc at z=0.05. 

In this section, we investigate the redshift dependence 
of the key morphological indicators in three bins of redshift 
within the MGC sample. Some care must be taken because, 
as we have seen, the disk fraction changes rapidly with lu- 
minosity; therefore, without first matching on luminosity 
we would have a higher disk fraction in the lowest redshift 
bins. We have broken the sample into low (0 < z <= 0.05), 
medium (0.05 < z <= 0.1) and high (0.1 < z < 0.15) redshift 
bins. Each galaxy in the 0.05 < z <= 0.1 bin was randomly 
matched to a galaxy in each of the other two redshift ranges 
with an absolute B magnitude within 0.03 magnitudes. The 
sample has 2169 galaxies in the 0.05 < ^ < 0.1 range. They 
were matched to a unique sample of 266 galaxies in the 
< z < 0.05 range and 1354 galaxies in the 0.1 < 2; < 0.15 
range. Each galaxy was weighted by the number of times it 
was matched. 

Figure [BT] shows the B/T distribution of the three dif- 
ferent redshift samples. These distributions are very similar 
in all redshift bins and our conclusions are unchanged if the 
field sample is taken as any of these bins. 

In ^3.3[ we claim that the observed difference in the 
mean asymmetries between the CN0C2 sample and the 
MGC sample is due to the better physical resolution of the 
CN0C2 images. To test this we reduce the resolution of the 
original image and the residual image of a representative 
sample of 60 galaxies from our CN0C2 sample by convolving 
with a Gaussian with different widths. Figure [B2l shows the 
asymmetry distribution of the original sample (thin, solid 
black line), and the asymmetries after broadening with a 1 
kpc (black, dashed line) and 2 kpc PSF (thick, solid black 
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Figure Bl. The B/T distribution of luminosity-matched samples 
of MGC galaxies in bins of redshift. The sample is divided into 
low (0 < z < 0.05; thin, dashed, black line), medium (0.05 < z 
< 0.1; thick red line) and high (0.1 < z < 0.15; thin, solid, blue 
line) redshift bins 



line). Clearly the asymmetry is reduced with poorer physi- 
cal resolution, which explains the higher asymmetries of the 
CN0C2 sample. For this reason, we do not compare the 
asymmetries between the two surveys and always use sam- 
ples matched in redshift within a given survey. 



APPENDIX C: IMAGES OF A SAMPLE OF 
GROUP AND FIELD GALAXIES AT Z=0.4 

In this section, we show a representative sample of the group 
and field galaxies from our CN0C2 z=0.4 sample. As dis- 
cussed in ^2.3[ the thumbnail images are from HST ACS 
observations. Each image is shown together with the gim2d 
model and the residual of the HST image after the model 
was removed. We show the group galaxies in Figure [CT] and 
the field galaxies in Figure IC2I In Figure IC3I we show im- 
ages of the nine group galaxies in the CN0C2 sample which 
are bulge dominated (B/T>0.5) and have high asymmetries 
(Rt + Ra >0.16). The first four of these galaxies show 
interaction features. 
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Figure CI. A representative sample of group galaxy images from the z=OA sample. For each of 24 galaxies is the HST ACS image (left 
panel), the GIm2d output model galaxy (middle panel), and the residual of the HST ACS image after the model is subtracted (right 
panel). Each galaxy is listed with the Bulge to total ratio (B/T) and Asymmetry (Asym) computed by GIm2d. 
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Figure C2. As in Figure [CTl but for a representative sample of 24 field galaxies in the z=OA sample. 
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Figure C3. The nine group galaxies in the z=OA sample which are bulge dominated (B/T>0.5) and highly asymmetric {Rt + Ra > 
0.16). The galaxies are shown with two stretches to show the interaction features. Prominent interaction features are seen in the first 
four of these galaxies. 
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Figure B2. The distribution of asymmetries (Rt+Ra) measured 
within 2 halflight radii for different physical resolutions for a sam- 
ple of 60 representative CNOC2 galaxies. The images were blurred 
with a Gaussian with a PSF of 1 (black, dashed line) and 2 kpc 
(thick, solid black line). The original, unblurred asymmetry is 
shown as the thin black line. 



